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Abstract. We provide a check of the accuracy of the auxiliary field formalism used to derive the Effective 
Hamiltonian for baryons in the Field Correlator Method. To this end we compare the solutions for the 
Effective Hamiltonian with those obtained from the solution of the spinless Salpeter equation. Comparing 
these results gives a first estimate of the systematic uncertainty due to the use of the auxiliary field 
formalism for baryons. 

PACS. 12.38.-t Quantum chromodynamics - 12.40.Yx Hadron mass models and calculations 



1 Introduction 



The advent of new ideas concerning quark-quark forces 
in QCD has led to a revival of interest in baryon spec- 
, troscopy. Various versions of the constituent quark model 
' [I] reproduce the octet and decuplet ground states but 
I have very different and even contradictory predictions on 
' the spectrum of excited states. It is therefore very im- 
portant to develop model independent methods that are 
. directly connected to the QCD Lagrangian and can help 
' in alternatively understanding baryon spectroscopy. 

One of such approaches is based on the Field Correla- 
tor Method (FCM) in QCD [g. PCM provides a promis- 
] ing formulation of the nonperturbative QCD that gives 
■ additional support for the quark model assumptions. The 
' application of this method for light mesons, heavy quarko- 
nia, heavy-light mesons and light and heavy baryons can 
be found in Refs. [3]. The key ingredient of the FCM is 
the use of the auxiliary fields (AF) initially introduced in 
order to get rid of the square roots appearing in the rel- 
ativistic Hamiltonian Q- Using the AF formalism allows 
one to write a simple local form of the Effective Hamilto- 
nian (EH) for the three quark system [6] , which comprises 
both confinement and relativistic effects and contains only 
universal parameters: the string tension a, the strong cou- 
pling constant as, and the bare (current) quark masses m^. 
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^ Historically the AF formalism was first introduced in [3] to 
treat the kinematics of the relativistic spinless particles. For 
a brief review of the AF formalism relevant to the problem 
considered in this paper see Sec. II of [5]. 



The EH has the form 



H 



3 

E 

i=l 



2^4 



(1) 



In Eq. ID), Ha is the non- relativistic kinetic energy oper- 
ator for masses /ii, V is the sum of the string potential 
VV(ri, 1*2 , Ta) and a Coulomb interaction term Vcouiomb 
arising from the one-gluon exchange. The string potential 
is 



Vy{yi, V2, rs) = (T) 



(2) 



where Tmin is the minimal string length corresponding to 
the Y-shaped configuration. Finally the /i^ are the oper- 
ator AF that have to be determined from the variational 
principle. 

Note that the the sum of the mass term and Hq in Eq. 
^ can be conveniently written as 
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2 fit 



f ) . (3) 



After taking the extremum of this expression in fii one 
ends with the standard relativistic kinetic energy operator 

i=l 

In this paper we use an approximate approach to con- 
sider the AF formalism first suggested in [7]. The AF are 
treated as c-number variational parameters. In this ap- 
proach one replaces the operators fj,i (r) depending on time 
parameter r by the c-numbers fii independent of r. The 
eigenvalue problem is solved for each set of nf, then one 
has to minimize (H) with respect to fii. Such an approach 
allows for a very transparent interpretation of AF: start- 
ing from bare quark masses rrii , we naturally arrive at the 
dynamical masses fii that appear due to the interaction 
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and can be treated as the dynamical masses of constituent 
quarks. 

An obvious disadvantage of the AF approach is that, 
as a variational method, it provides only an upper bound 
to the mass spectrum. So far the accuracy of this approx- 
imate solution for relativistic systems has been checked 
numerically only for mesons [HH]. The principle objective 
of this work is to test the AF method for baryons. We 
implement the AF method to calculate the baryon masses 
and then perform similar calculations using the relativistic 
Hamiltonian 



H = 



V. 



(4) 



Although being formally simpler the Hamiltonian ([T]) is 
equivalent to ^ up to the elimination of the AF (see e.g. 
Ref. [S]). We refer to an eigenvalue equation with Hamilto- 
nian dl]) as the spinless Salpeter equation (SSE). In QCD, 
it arises from the Bethe-Salpeter equation replacing the 
interaction by the instantaneous potential V and consid- 
ering a limited Fock space containing qqq states only. 

In this paper, we study the confinement plus Coulomb 
energies for the ground 5-wave and orbitally excited P- 
wave states of nnn, nns and ssn baryons and disregard 
the spin dependent forces, which are not relevant for our 
consideration. 

The baryon masses in the AF approach are calculated 
using the hyperspherical method, while those in the SSE 
are calculated variationally. The numerical algorithm to 
solve the three-body problem variationally is based on an 
expansion of the wave function in terms of harmonic oscil- 
lator functions with different sizes [10] . The details of tech- 
nical aspects can be found elsewhere [11]. It was proved 
to give results of good accuracy if the expansion is pushed 
sufRciently far (let say up to 16-20 quanta). Moreover it 
can deal easily either with a non-relativistic or relativistic 
expression for the kinetic energy operator. 

We find an accuracy of the AF method for hyperons to 
be about 6 % at worst, which is quite reasonable to justify 
application of the AF formalism. 

The paper is organized as follows. In Sec. 2, we briefly 
review the EH method. The application of this method 
for the baryons was described in detail elsewhere [l2l[T3]. 
Here we give only a brief summary important for our par- 
ticular calculation. In Sec. 3, we discuss the hyperspherical 
approach, which is a very effective numerical tool to solve 
this Hamiltonian. In Sec. 4, we provide a few numerical 
examples illustrating the accuracy of the hyperspherical 
solutions. In Sec. 5, predictions of the AF method are 
compared with those obtained from the solution of the 
spinless Salpeter equation (SSE). Section 6 contains our 
conclusions. 



^ Here and below the symbol n stands for the light quarks u 
or d. 



2 The baryon masses in the AF method and 
SSE 

The baryon mass in the FCM is given by 

M^^ = Afo^^ + C^^, (5) 



Mr 



3/2 

V ( ^ -I- — 



i?o(/i.) (6) 



where EQ{fii) is an eigenvalue of the Shrodinger operator 
Hq + V, the constant AF fii are defined from the minimum 
condition 

9Mo^^(m„M.) 



9 Ml 



= 0, 



(7) 



and C^^ is the quark self-energy correction which is cre- 
ated by the color magnetic moment of a quark propagating 
through the vacuum background field [14] . This correc- 
tion, which can be added perturbatively, adds an overall 
negative constant to the hadron masses: 

C^^ = --E^' (8) 



where 1/Tg is the gluonic correlation length. In what fol- 
lows we use Tg = 1 GeV. 

The function 77 (t) is defined as 

/•OO 

r]{t) = t / Ki{tz)e-'' dz, (9) 
Jo 

where Ki is the McDonald function. A straightforward 
calculation yields [2] 



1 + 2t^ 

(1-^2)2 ~ (l-t2)5/2 

1 + 2<2 3^2 

(1-^2)2 " (t2 _ 1)5/2 



In , t < 1, 



arctan(Vi^ - 1), t -^m) 



Note that 77(0) = 1 and rjit) ~ 2/<2 as t ^ 00. 

The baryon mass in the SSE approach is given by 



M 



SSE 



= M, 



SSE 



C 



SSE 



(11) 



where Mq^^ is an eigenvalue of the relativistic Hamilto- 
nian ([U) and the C^^^ are given by ([5]) with the obvious 
substitution fii — > Ui, where 



(12) 



are the average kinetic energies of the current quarks. 

We will not perform a systematic study in order to 
determine the best set of parameters to fit the baryon 
spectra. Instead, in what follows we employ some typical 
values of the string tension a and the strong coupling con- 
stant as, which have been used for the description of the 
ground state baryons [T^]: a = 0.15 GeV2 and = 0.39. 
In our calculations we use the values of the current light 
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quark masses, rUu — md — 9 MeV, and rus — 175 MeV. 
As in Ref. [12] we neglect the spin dependent potentials 
responsible for the fine and hyperfine splittings of baryon 
states. 

Our aim is to compare the baryon masses given by Eqs. 
([5]) and (fTT|) . To this end we first solve the non-relativistic 
Schrodinger equation with the confining and Coulomb in- 
teractions to determine the constituent quark masses 
and the baryon masses . Efficient methods to deal 

with the Y-shape interaction rely either on Monte-Carlo 
algorithms ^15,16^ or the hyperspherical method 17J. We 
use the latter approach. 



3 Outline of the hyperspherical formalism. 

In this section, we briefly review the hyperspherical method, 
which we use to calculate the masses of the ground and 
excited hyperon states. 

The baryon wave function depends on the three-body 
Jacobi coordinates 




(13) 



(i,j,k cyclic), where and /iy^fc are the appropriate re- 
duced masses: 



My 



i-^ij, k 



Ml 



(14) 



and Mo is an arbitrary parameter with the dimension of 
mass, which drops out in the final expressions. There are 
three equivalent ways of introducing the Jacobi coordi- 
nates, which are related to each other by linear transfor- 
mations with the Jacobian equal to unity. In what follows 
we omit the indices i and j. 

In terms of the Jacobi coordinates the kinetic energy 
operator _ffo in HI) is written as 



1 



1 b_d_ 
2^ \dB^ ^ RdR 



(15) 



where R is the six-dimensional hyperradius that is invari- 
ant under quark permutations, 

p^RsinO, \ = Rcos9, < 6* < 7r/2, (16) 

n denotes five residuary angular coordinates, and i^(]7) 
is an angular operator 



1^ 



whose eigenfunctions (the hyperspherical harmonics) sat- 
isfy 

L'{n)Y[K]{e,np,nx) = -K{K + A)Y[K]{B,np,nx), (18) 

with K being the grand orbital momentum. 

The wave function "(/"(p, A) is written in a symbolical 
shorthand as 



V'(p,A) = ^V[K](i?)V](^). 

[K] 



(19) 



where the set [K] is defined by the orbital momentum of 
the state and the symmetry properties. 

We truncate this set using the approximation K = 
Kmin- We comment on the accuracy of this approximation 
latter on. Our task is then extremely simple in principle: 
we have to choose a zero-order wave function correspond- 
ing to the minimal K for a given L (-fCmin ~ for L — 
and Kjriin — 1 for L — 1). The corresponding hyper- 
spherical harmonics are 



Yo = 
Y„ = 



, K = 0, 



3 



6 A 



7r3 R 



K = I. (20) 



For nns baryons we use the basis in which the strange 
quark is singled out as quark 3 but in which the non- 
strange quarks are still antisymmetrized. In the same way, 
for the ssn baryon we use the basis in which the non 
strange quark is singled out as quark 3. The nns basis 
states diagonalize the confinement problem with eigen- 
functions that correspond to separate excitations of the 
non-strange and strange quarks {p- and A excitations, 
respectively). In particular, excitation of the A variable 
unlike excitation in p involves the excitation of the "odd" 
quark (s for nns or n for ssn). The nonsymmetrized uds 
and ssq bases usually provide a much simplified picture 
of the states. The physical P-wave states are neither pure 
SU(3) states nor pure p or A excitations but linear combi- 
nations of all states with a given J. Most physical states 
are, however, closer to pure p or A states than to pure 
SU(3) states [H]. Note that for the nnn baryons, the p 
and A excitation energies are degenerate. 

Introducing the reduced function u^{R) 



(21) 



where 7 = for L ~ 0, 7 = p, A for i = 1 0, the new 
variable 



X = y/jM)R 



\ " Ml M2 2 



1/2 



M 



M2 M3 J2 

M 



' 23 



M3 Ml ^2 

M 



(22) 



''' In what foliows, for ease of notation we wili drop the mag- 

(17) 

netic quantum numbers of the vector spherical harmonics. 
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and averaging the interaction V = Vy + Vc over the six- 
dimensional sphere i7 with the weight one obtains 
the one-dimensional Schrodinger equation for Uj{x) 



2 ( Eo 



— Vjyu., J u.^ 



V.,{x) u^{x) = 0(23) 



where V^(x) = [x) + V^^^y^^^{x), 



V^{x)^ I \Y^{9,x)\^VY{r,,r2,r3)dn = 



a hi, 



(24) 



and 



^Coulomb (^) = - o"^ 



2 a. 



J - . ■ fij 



-as — VMo- 

3 X 



In what follows we denote 



(25) 



(26) 



Then the straightforward analytical calculation of the in- 
tegrals in ((25)) yields 



16 



ao 



^-3^1^+VrT^j^' (27) 
32 



15 TT 

32 



V2 



K 5k + 6 

1 + K 1 + K 



K 4 + 5k 



l^TT \y/2 y 1 + K 1 + K 



VJ^, (28) 



Vm- (29) 



For K = 1 (the nnn system) Op = ax- The correspond- 
ing expressions for bj are more complicated (see, e.g., the 
appendix of Ref. jlj). 



4 Accuracy of the hyperspherical 
approximation 

A few words concerning the accuracy of the approxima- 
tion K — if mill are in order. An illustration of the accuracy 
of the hyperspherical approximation K = i^min is given 
by the results presented in Table [T] This Table compares 
the eigenvalues Eq in Eq. ^ for the nnn, nns and ssn 
systems obtained using the variational method and those 
calculated from Eq. ([23]) with K = Kmin El- In all cases 
the dynamical masses fii are the same as were found from 
the minimum condition ([7]) for the Y-shaped string poten- 
tial [13] . For technical reasons the variational calculations 



Recall that, as was stated in Sec. 3, the p and A excitation 
energies for the nnn baryon are degenerate. 



have been performed not for the genuine string junction 
potential but for its approximation by a sum of the one- 
and two-body confining potentials [19j 



1 



[Va + Vcm), 



(30) 



where Va with the sum of the two-body confining poten- 
tials is 



(31) 



and Vc is the sum of one-body center-of-mass string po- 
tentials: 



i 



(32) 



{i,j,k cyclic), where Rem is the center-of-mass coordi- 
nate. Table [1] also compares eigenvalues E^ for the gen- 
uine string potential Vy with those for the confining po- 
tentials Va, Vqm and Vm with the same string tension. 
The confining potential Vcm overestimates the eigenval- 
ues of the genuine string junction E^ while the potential 
Va underestimates the E^ , i.e. E^ < E^ < E^^^ (com- 
pare columns 6, 7 and 8 of Tabled]). The values of the two 
columns 7 and 8 are in reasonable agreement with the ref- 
erence results of column 6. In line with expectations jl9j . 
the eigenvalues for the genuine string junction change lit- 
tle if we use Vm instead of Vy. Simulation of the genuine 
string junction potential by a sum of the two-body con- 
fining potentials ([50]) (column 9) is a good approximation 
in all cases: using Vm results in a ~ 20 MeV or 1 — 2 % 
downwards shift of Eq for all states (compare columns 6 
and 9). Let us note that Hamiltonian[T]with potential Vm 
gives eigenvalues wich are, to some MeV, the arithmetic 
mean of the eigenvalues with potential Va and Vcm- So 
the contributions of Va and Vcm to Vm are nearly evenly 
distributed. 

The last column 10 contains the eigenvalues E^^^^^ 
calculated using the variational method briefly described 
in Sect. 1. Comparing the column 9 and 10 of Table [1] we 
conclude that the hyperspherical and variational results 
are close enough to validate the approximation K = Kmin- 



5 Comparison of the AF and SSE results 

Table [2] compares the baryon masses computed using the 
AF and SSE formalisms. In this Table we list the masses 
of the nnn, nns and ssn states with L — 0,1. The en- 
tries labeled AF have been calculated from Eq. ([23]) with 
K ~ -fsTmin, while the entries labeled SSE have been cal- 
culated using the variational method for the relativistic 
Hamiltonian (|4]). In both cases, we approximate the Y 
shaped string potential by the expression ([50)) . As was 
mentioned in the Introduction the comparison of the AF 
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results with those evaluated from the solution of SSE has 
been performed only for the qq mesons with the conclusion 
that the variational AF method gives a systematic over- 
estimation of order 5-7 % [lUH]. Our calculations show the 
similar results: the relative deviation 



e = 



Mr 



MSSE 



(33) 



is positive and for most considered states does not ex- 
ceed 6% 0. The accuracy of the AF approach does not 
seem to be very sensitive to the bare light-quark masses. 
The quantum numbers of states have a stronger influ- 
ence on the accuracy. In particular, e for the L — 1 
states are uniformly smaller than those for the L = 
states. Curiously, the self-energy corrections C^^^'-^ and 
j[fSSE g^gj.gg even with better accuracy (typically within 
5% or even better) in spite of the fact that the differ- 
ence and uji in some cases (e.g. for the A excitation in 
the ssn) comprises 30%. As for the excitation energies, 
A = Mb{L = 1) - Mb{L = 0) evaluated using the 
AF and SSE methods, they practically coincide for the 
ssn baryons and differ no more than 30 MeV for the 
nns baryons. Taking into consideration that we neglect 
the spin interactions the baryon energies calculated us- 
ing SSE agree reasonably with the data [20]. For instance, 
for L = we get 1{N + Zi)theory = 1062 MeV versus 
i (TV + Zi)exp = 1085 MeV and i(yl + ^7 + 2 i:*)thcory = 
1220 MeV versus 1{A + S + 2E*%^p = 1267 MeV. A 
similar correspondence exists for the other states consid- 
ered in this work. 



6 Conclusions 

In this paper we have tested the quality of our previous 
study of the masses of the S- and P- baryon states ob- 
tained in the FCM with the use of the AF formalism. To 
this end we have compared the AF results with those ob- 
tained from the solution of the SSE with the same interac- 
tion. The main purpose was to check whether the results 
obtained within these two methods are similar. We have 
found that they agree within ~ 100 MeV for the absolute 
values of masses and with much better accuracy for the 
excitation energies. Thereby our study supports the AF 
basic assumptions by the compatibility of its mass pre- 
dictions with the masses derived from the SSE. Moreover, 
this comparative study gives better insight into the quark 
model results, where the constituent masses encode the 
QCD dynamics. 
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Table 1. Comparison of the eigenvalues Eo of the Hamiltonian Ho + V in Eq. ^ for the baryon ground states and the p 
and A excitations obtained from the hyperspherical solution of Eq. (|23|l {E^) and variational solution (-Eo^tmr)- See the text for 
further explanation. 



Baryon 


L 


Excitation 


Ml 






Eo 


Eo 


Eq 


-^0 var 


nnn 



1 


p, A 


408 
457 


408 
457 


1318 
1638 


1366 
1697 


1230 
1532 


1299 
1615 


1297 
1612 


nns 



1 
1 


P 
A 


414 
482 
441 


453 
459 
534 


1291 
1611 
1614 


1339 
1670 
1676 


1204 
1506 
1508 


1272 
1589 
1593 


1271 
1587 
1591 


ssn 



1 
1 


P 
A 


458 
520 
483 


419 
424 
506 


1266 
1592 
1588 


1313 
1653 
1646 


1181 
1487 
1485 


1248 
1571 
1567 


1248 
1569 
1566 



Table 2. Comparison of baryon masses calculated using the AF approach and SSE. The symbol Vi denotes either the constituent 
quark masses /ii or the average kinetic energies of the current quarks lOi. Shown are the masses Mq and M§ without the 
self-energy corrections, the self-energy corrections C^^ and C'^^^ , M = Mo + C (all in units of MeV), and the relative error 
£ defined by Eq. ((33)) . 



Baryon L Excitation Method vi ^ V2 Mo C M e(%) 



nnn 







AF 


408 


408 


1911 


-702 


1209 


13.8 








SSE 


394 


394 


1788 


-726 


1062 




nns 







AF 


414 


453 


1946 


-648 


1298 


6.4 








SSE 


396 


484 


1877 


-657 


1220 




ssn 







AF 


458 


419 


1182 


-598 


1384 


6.3 








SSE 


404 


465 


1904 


-602 


1302 




nnn 


1 


p, A 


AF 


457 


457 


2301 


-627 


1674 


9.1 






p, A 


SSE 


440 


440 


2186 


-651 


1534 




nns 


1 


P 


AF 


482 


459 


2356 


-581 


1751 


6.0 






P 


SSE 


464 


465 


2245 


-594 


1652 








A 


AF 


441 


534 


2330 


-592 


1738 


1.5 






A 


SSE 


415 


592 


2315 


-603 


1712 




ssn 


1 


P 


AF 


520 


424 


2362 


-552 


1810 


4.1 






P 


SSE 


478 


530 


2302 


-564 


1738 








A 


AF 


483 


506 


2367 


-540 


1827 


4.4 






A 


SSE 


503 


391 


2295 


-545 


1750 





